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The theory of ortho to para hydrogen molecule catalysis on 
magnetic surfaces was carefully reviewed during the course of this 
study. The theory was developed from 'first principles’ by deriving a 
master equation, which was then evaluated in a form which includes 
the effects of the closeness of approach of the H2 molecule to the 

surface, the physical state of the molecule on the surface (gas or 
liquid), and the magnetic state of the catalytic surface: para- 

magnetic or having collective magnetic excitations due to long range 
magnetic interactions (ferromagnetic or antiferromagnetic). 

The theory was then brought to a numerical stage and evaluated 
for a variety of models, which were reported in six short reports 
which were sent to NASA/Lewis. Copies of each of these reports are 
included as the body of this final report. These reports contain 
numerical evaluations of the absolute rate of para to ortho hydrogen 
catalysis for a variety of dimensionless parameters which are 
useful for describing the rate of catalysis for different materials, 
and different catalytic conditions. See report VI for the best 
evaluation ana testing of a numerical model (using transition state 
theory for the migration of hydrogen molecules between catalytic 
sites) written in VAX FORTRAN 77. 

The final conclusions of this study are briefly summarized as 
follows: regardless of the detailed structure of any catalytic 

material, the rate of para to ortho hydrogen on a (magnetic) 
catalytic surface is primarily a function of only two parameters: 

(1) The closeness of approach of the H2 molecule to the 

magnetic moments on the catalytic surface: the catalytic 

rate is a maximum for a minimum distance of approach of the 
hydrogen molecule to the surface layer containing the magnetic 
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moments. A surface layer containing 'bare magnetic moments' 
would be the best in regards to this criterion. 

(2) The surface lifetime. The surface lifetime must be as short 

as possible, or as close as possible to lO' 1 ^ sec. as can be 
physically achieved. This lifetime parameter (essentially the 
surface lifetime) is difficult to control and may be nearly 
impossible to measure. It aopears to be difficult to achieve so short 
a surface lifetime at the target temperature of 120 K for the 
hydrogen gaseous phase. Studies of the para --> ortho catalytic rate 
versus lifetime showed considerable degradation of catalytic 
performance as the surface iifetime rises from its optimum value. 

The major differences between different catalytic surfaces were 
deemed, in the conclusion of this study, to be related to the physical 
nature of the surface as regards its molecular surface complex- 
ity and the extent of physical binding of the H 2 molecules to 
the surface. 

The molecular surface complexity may be considered to be those 
features of the surface which prevent optimum approach of the 
hydrogen molecules to the active catalytic sites. The lifetime of 
the molecule on the surface is related to the nature of the peaks 
and valleys of the adsorption potential that a particular molecule 
would ’see’ as it migrates (by surface diffusion) from one peak of 
suface potential energy (surface atom or molecule) to another. So 
the only way a really detailed understanding of this surface 
catalytic process can be achieved is tc model the physical potentials 
which weakly bind the H 2 molecule to the surface. This requires 

structural knowledge of the nature of the surface of the catalytic 
material. Building such information into a model for surface 
diffusion and binding would then give a more complete picture of the 
overall process of para to ortho hydrogen catalysis, which could then 
be used to correlate the relative catalytic activities of different 
materials. Some experimental knowledge of the magnetic 
character of the catalytist is also considered important for 
empirical correlation of catalytic activities for different materials. 
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ABSTRACT 


The theory of Petzinger and Scalapino (Physical Review B, 8, 2 66 
(1973) ) has been thoroughly reviewed, and all of the basic equations 
for paramagnetic para to ortho hydrogen catalysis re-derived. The 
paper has only a few minor phase errors and errors of omission in 
description of the theory. Three models (described by Petzinger and 
Scalapino) for the rate of para to ortho h’ 2 catalysis were worked out. 
and uniform agreement obtained to within a constant factor of 'In’. 
The analytical methods developed in the course of this study were then 
extended to two new models, which more adequately describe the process 

of surface catalysis including transfer of hydrogen molecules onto and 
off of the surface. 


All five equations for the para to ortho catalytic rate of 
conversion are described here. The two new equations describe the 
catalytic rate for these models: 

H 2 on the surface is a two-dimensional gas with lifetime 

H 2 on the surface is a two-dimensional liquid undergoing Brownian 
motion, (diffusion) with surface lifetime 'x'. 

A numerical study of these rate equations will be the subject of the 
next report. 


Contents : 

I . General Methodology 

II. Models 

III. Analytical Results for Each Model. 

IV. The Next Step: Numerical Studies. 
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I . General Methodology 

There are two papers in the literature which appear to have a 
reasonable physical formulation for the description of both the 
paramagnetic nature of para to ortho hydrogen catalysis and the 
characteristics of the solid state surface of the catalyst. 

These are: 

K. G. Petzinger and D. J. Scalapi.no, 'Para to Ortho Hydrogen 
Conversion on Magnetic Surfaces', Phys . Uev. 8, 266 (1973). 

Y. Ishii, 'Theory of Non-Dissociative Ortho-Para Conversion on 
Magnetic Surfaces', Progress in Surface Science, 21, 163-203 (1986). 

These papers present an overall view of the theory of this 
catalysis process, and some of their results. However, a great deal of 
the methodology is omitted (essentially no appendices), sc the only way 
to acquire the technological ability to APPLY their results is to 
retrace their 3teps and learn the physical-mathematical methods which 
led to their conclusions. A good part of this process has now been 
accomplished, which has led to some new results which may be of use in 
the practical matter of designing good catalysts for para to ortho 
hydrogen conversion. 

The general description adopted in this study concerns only the 
lowest para state and lowest ortho states(s) . There is only one lowest 
para state, but the lowest ortho state has nine components (degeneracy 
- 9 ), since for this state I-i and L-l. (3X3 - 9 ) 
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The average frequency of the transition (in radians/sec) is designated 
by the symbol u Qp . So the energy of the transition is 

E(ortho) - E(para) » ttu 

op 

The practical results are found by integration and summation of 
terms xn a master equation which describes anv possible catalytic 
system of hydrogen molecules of either para or ortho type. Derivation 
and proof of the correctness of the many terms in this master equation 
has. up to now. been a major object of this study. 

The first step is to derive, for any model, a relaxation function 
designated by fl)(t), which is a function of time: 


Each unique model < > a unique <0(t) 





The rete of catalysis is then found by integrating the relaxation 
function over all time: 
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The Fourier Transform coefficients. 0^ (l) (q). have been calculated and 
agreement obtained (except for some phases) with both earlier authors 

quoted above. These terms are critical for the purpose of developing a 
reliable theory. 
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For each model, then, one must evaluate: 


The orientation and selection rule term: 

A 

(b - orientation of the* molecule) . 
(integrated over all orientations) 


C 


The correlation function for motion of the hydrogen molecules across 
the surface: 
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The correlation function for the motion of the electronic spins on the 
surface: 



CO 

(o) 



b |ec . 


The electronic spins are the actual catalyzing sites for the conversion 
of B 2 nuclear spins from para (paired) to ortho (unpaired). In the 
calculations described here, this term only measures the magnitude of 
the electron spins and their concentration (paramagnetic case) . 

T.n addition, one must evaluate (once and for all) the Fourier 
Transforms of the dipolar form factors, which give the coupling between 
the dipolar interaction as a function of spatial coordinates and the 
extended wave mo ; ons of the lattice described as function of che wave 
vectors q. Only h f of the first Brillouin zone is considered, (q « 0 
to 2»/oi where a is the lattice constant of the surface) . 


Once each of the terms above has been evaluated for each model, 
then the catalytic rate is obtained by an integration over time and a 
summation over the wave vector q. A standard method of solid state 
physics is then used to convert the sum over q to an integral. Thus, 
all expressions presented here as analytical equations are essentially 
integrals over the q wave vector, which must be determined numerically. 

II. The Models. 

Each model described here is for catalysis on a paramagnetic 
surface, with a certain concentration (N s /N) of paramagnetic sites in 
the surface. 

Model PS-l: The H 2 exists on the surface as a two-dimensional gas 
with infinite lifetime. (H 2 molecules never exit from the system). In 
time interval 0 to t, such a molecule moves a distance over the surface 
'r', where r is found by knowing the RMS speed of a molecule derived 
from a 2-diroension Boltzmann distribution of molecular speeds. 

Thus : 


lr(t)-r(0)|«v t 

rms 


Model PS-2: The H 2 exists on the surface as a two-dimensional 

liquid with the H 2 molecules executing a Brownian motion random walk 
described by a diffusion coefficient 'D' . Thus, by the usual relations 
for the mean squared displacement of a diffusing particle (in a two 
dimensional space) the displacement after time t i.s: 
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lr(t) - r (0) I - SQRT( 4Dt ) 

Thus there is a different dimensionality in the integrand of the final 
equations to be integrated, between these different models. 

Model PS-3: The H 2 exists, above the catalytic surface, as a 

dilute gas and simply hops on and off of the surface. The catalysis 
only occurs on the surface, during the lifetime of a hydrogen molecule 
on the surface of -c seconds. 


T h esq HfiW model s* worked out new with the methods used to derive the 
first three above, are more likely to describe the real-world para to 
ortho hydrogen catalysis situation for paramagnetic catalvsis: 

Model I: The H 2 molecule exists as a two-dimensional gas on the 
surface, and hops on and off of the surface into an external extended 
gaseous phase with temperature T. The lifetime of an H 2 on the surface 
is r seconds. 

Model II. The H 2 molecule exists on the surface as a two 
dimensional liquid, and hops on and off of the surface into an external 

gaseous phase with temperature T. The lifetime on the surface is again 
t seconds . 



V 


The process of model building then consists of converting 
simple picture, given above, into a mathematical description. 

The final step in the process is to express the rates of para to 
ortho conversion in terms of dimensionless parameters which may then be 
used to characterize different catalyzing surfaces. The magnetic and 
structural differences in catalytic rates between different surfaces 
may then be expressed in terms of the experimental values of these 
dimensionless parameters. 
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III. Analytical Results for Each Model 


Model PS-l: h 2 on surface is a Two-dimensional gas 

(Surface lifetime t = infinite) 


<H 2 in gaseous phase) 


Rate - zero 
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Model PS-2: h 2 on. the Surface Exists as a Two-dimensional Liquid 

(Surface lifetime t is infinite) . 


( H 2 in gaseous phase) 

4 I 
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Model PS-3: 'Jump' Model. (Surface Lifetime - t seconds) 
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( in gaseous phase ) 
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where VT 

p-~ >o 


rate of para to ortho conversion (sec -1 ) 


Parameter 9 - M op T ls dimensionless . 

Rate depends parametrically on 0. 

There is also some dependence on the upper limit of the integral: 2q c z, 
but the main effect is expressed through the dimensionless parameter 9. 
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Practical Model 


I: H 2 on the Surface exists as a two-dimensional 
gas with finite surface lifetime m x seconds . 


H, 


(Gaseous Phase) 
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Practical Model II ; 


Hj exists on surface as a two-dimensional liquid 
with surface lifetime - t seconds. 


(Gaseous Phase) 

k i 


Rate ~ l/r 
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Practical Model II. (Continued). 
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a o " ^ att i c e constant for the surface lattice. 

Z - the height of the H2 molecule above the surface in the 
adsorbed two-dimensional layer* 
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Definition of the universal catalytic constant K: 

HL f \ limj ( b ') mp pn ( 1 \ 

t InJ 3 \ Qa I \ 28 / 
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Ns/N 


P, 


2B 


the fraction of surface atoms which carry 
unpaired electronic spins. 


'Spin' of the paramagnetic atom on the surface. 


e.g. 1/2, 1,3/2.. 


bond length of the H 2 r. olecule. 


lattice constant of the surface layer 


‘p* "A “ magnetic moments of the proton or the paramagnetic 
atom . 


the height of the H2 molecule above the surface plane 
containing the paramagnetic catalyzing atoms. 


the rotational energy constant of the H2 molecule 
when present in the two-dimensional surface layer 
(may be different from the gaseous phase) . 
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IV. The Next Step: Numerical Studies. 

The next step is to generate curves, by numerical integration, of 
the catalytic rate versus values of the dimensionless parameters 
characteristic of each model. The comparisons will be: 

1. Jump model versus 6. 

2. Model I versus dimensionless parameters o and e. 

3 . Model II versus dimensionless parameters r and p . 

The general expectation is that the families of curves so generated 
will, in each case, show a peak in catalytic rate for some values of 
the parameters. Thus, in some cases an increase in a parameter value 
may lead to an increase in catalytic rate, and in other cases an 
increase in a parameter value may decrease the catalytic rate. So it 
will be very to determine which side of the peak is a given catalyst. 

The numerical comparison of the models above may yield some more 
quantitative information as to optimal types of physical systems. 



The NEXT STEP: 

The Ferromagnetic and Antiferromagnetic problems will be the next 
step. These equations will probably be similar to the ones above, 
except for resonance terms representing resonant exchange of energy 
with the many body magnetic excitations, or magnons . Tnis may cause 
greater catalytic rates than for the paramagnetic case, in certain 
cases . 
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Theory of Ortho to Para Hydrogen Catalysis 
II. Analysis of Three Models Due to Petzinger and Scalapino. 

August 10. 1989 

R. E. Coffman 
Chemistry Department 
University of Iowa. 

ABSTRACT 

Three simple models for para to ortho hydrogen catalysis on 
paramagnetic surfaces, originally due to Petzinger and Scalapino. have 
been numerically simulated. These models are: A two-dimensional 

liquid on the surface, a two-dimensional gas on the surface (each with 
infinite surface lifetime) , and a dilute gas which hops on and off of 
the surface with lifetime *t'. The integrals determining the 

functional dependence on certain macroscopic dimensionless parameters 
have been evaluated with the IMSL numerical routine QDAGS (Gauss-Konrad 
rule), and then plotted versus the dimensionless parameter character- 
istic of the rate of catalysis for each model. The results are 
presented in a series of graphs for several ranges of the dimensionless 
parameter for each model, and for different values of the height of the 
H 2 molecule above the surface. One cannot generally say that one model 
necessarily leads to higher catalytic rates than another, but the two 
dimensional gas model clearly peaks at smaller parameter values than 
the other models and attenuates faster versus its dimensionless 
parameter. Thus, the two-dimensional gas model leads to larger para to 
ortho catalytic rates over a smaller parameter range than the other two 
models considered. 


CONTENTS 


I • Analytic Methods 

II. The Models. 

III. Intermediate Equations for the Rate of Para to Ortho 
Hydrogen Conversion. 

IV. Final Results: Condensed Equations for rate of Catal 

V. Interpretation and Conclusions. 
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VII. Plots. 


I. Analytic Methods. 




The rate of para to ortho catalysis of hydrogen molecules in 
contact with a paramagnetic surface is given as the integral of a 
relaxation function: 
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Evaluation of this functions requires evaluation of the orientations 
terms, the correlation functions and the Fourier Transform coefficients 

here indicated as ^(q). A.U of these steps have been independently 
performed and checked for accuracy. 



II. Models and Assumptions. 


It has been established that these assumptions are correct 
(within a certain degree of approximation) for the problem at hand: 

(A) The orientation terms are averaged over all orientations of the H 2 
molecules on the surface. The result is: 


< G^b) > = 4n/3 


(B) The correlation functions for the motion of the B 2 molecule on and 

off of the surface were evaluated for three models. The results for 
each model are as follows: 
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2 Dim-gas on surface 
* 3 infinite 


-Dq 2 t 


2 Dim-liquid on surface, 
* 3 infinite 


~t/x 


Dilute gas on surface. Jump 


Model . 


* 3 finite. 
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(C) The correlation function for the motion of the electronic spins 
on the surface is evaluated for the case of non-interacting spins, 
with surface concentration N /n. The result is: 

W 


cn c<> 
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(D) The sum over all of the Fourier Transformed 'beta' coefficients 
is used. This sum is: 
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which follows from the detailed evaluation of the Fourier Transforms of 


the dipolar form factors expanded in the lattice wave vectors a. 
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III. Intermediate Equations for the Rate of Para to Ortho Conversion. 
A set of equations, which involve double integrations (after 

-k — i 

transforming the sum over q into an integral over q) , are found for 
each of the three models considered here. The intermediate equations 
which must be evaluated in order to find the total rates of para to 
ortho conversion are as follows: 
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IV. Final 
catalysis . 

MODEL PS-i : 
where 

and 

MODEL PS-i- 


Results: Condensed equations for the rate of para to ortho 


H 2 “ 2 Dimensional gas on the surface, 
residence time z - infinite. 



KF ( a ) 


F ( a) 



a 



2 - 




( dimensionless ) 


H 2 " 2 Dimensional liquid on the Surface, 
residence time r * 


infinite . 



MODEL PS 


where 


and 


4E> 


( dimensionless ) 


H 2 * dilute gas; Jump model. 

Molecules jump on and off the surface, 

with residence time on the surface = r seconds. 
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The common coefficient ' K ’ is, for each case: 


+tr f Hi\ 
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which has a strong dependence on Z, i.e. K - l/Z 6 . 


Also: 


Ng/N - fraction of surface sites (N> occupied by paramagnetic 

atoms (N > . 
s 

S - the spin of a surface paramagnetic atom. 


The integral limit is approximated: 


q_Z ■ n( z/a ) 
c o 


which follows from: 


q c ~ n/a Q - cut-off value of the wave-vector; limit of first 

Brillouin zone (i.e. the shortest possible 
wavelength equals the lattice constant) . 


and the assumption that the distance of the H 2 from the surface is 
approximately equal to the surface latticf constant. a Q . Tnus: 

Model PS-l: Integral limit - 2c 3 c z ■ 2n(Z/a 0 > 


Model PS-2: 


Model PS-3: 


- q c Z/2 - 0 . 5n ( Z/a Q ) 

- 2q Z - 2 n ( Z / a ) 

C O 


These limits exert considerable influence on the values of the 
integrals which control the rate of catalysis for each respective 


model . 
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V. 


Interpretation and Conclusions . 



The integrals F(a), G(y) and H(6) were numerically evaluated using 
the IMSL program QDAGS * which is appropriate for an integral with 
finite limits. (The interval IA.B] is subdivided into subintervals; 
each subinterval is estimated with a 21-point Gauss-Konrad quadrature 
rule). The accuracy was set to about six figures (single precision; no 
difference was noted on switching to double precision) . 

The FORTRAN code for each integral was written into a single 
program which then put the data for the three functions versus the 
dimensionless parameter (for each) into a single file for plotting. 
The accompanying figures show, on the same plot, a numerical comparison 

Of: 


F(o) 

versus 

a . 

i 

G(y) 

versus 

r . 

j 

H(0) 

versus 

0 . 

i 

i 


•I 

The first six figures: 

G-l through G-6 

are plots versus a (or y or 9) over the range 0 to 3 , for different 
heights 'Z* of the H 2 molecule above the paramagnetic surface, measured 
in units of the surface lattice constant a Q . The range of values is: 

Z/a Q “ 0 . 2 , 0.5, 1 . 0 , 2 . 0 , 5 . 0 , 10.0 
(Figure G-3 is the same as Figure 2 in Petzinger and Scalapino's pa.'er, 
with, evidently, good numerical agreement!) 



The following three figures: 


H-i through H-3 

are plots versus the independent parameter a (etc.) over the range 0 to 
10 and for heights 'Z' having values: 

Z/a 0 - 0.5, l.o, 2.0 

The last three figures: 

I-i through 1-3 

are plotted over the range of the independent parameters a (etc.) over 
the range 0 to 30. and for the same range of Z/a Q as for the previous 
three figures . The purpose of this last set of plots is to emphasize 

the peaked nature of the catalytic rate for the two-dimensional gas 
model . 

Since all results are derived from the same master equation 
(accurate to all factors of 2. * , 'e', and the Uke) it seems 

reasonable to draw some limited conclusions already from the comparison 
of these three models for the different parameter ranges and heights of 
the Hj layer above the catalytic surface. 

» 

aome Conclusions: 

The two-dimensional gas model peaks over a smaller range of its 
independent parameter than the other two models. This property is 
traceable to the presence of an exp(-a 2 ...) terms in the integrand for 
this model, which does not occur in such a manner (not as an 
exponential!) attenuating the rate controlling integral for the other 
two models. (See figures 1-1 through 1-3). 
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Another conclusion is that a practical model for para -> ortho 
hydrogen catalysis must incorporate the rate of transfer of the H 2 
molecules into the over-all description of the catalytic rate. Such a 
description of the rate of catalysis would then be expected peak when 
plotted versus the parameter a or r» as well as a peak when plotted ; 

versus another dimensionless parameter containing the inverse of the 
surface lifetime. 1 /t. 


The formal difference between the two dimensionless parameters 
used for calculating the para to ortho hydrogen conversion rates for 
the first two models: 


and 


UJop 5* 


'll, 


n»nS 


*2" 

Tx” 


n 

also accounts for some of the differences between the shapes of the 
rate curves for the two-dimensional gas case compared to the case of a 
two dimensional liquid, each with infinite residence time. (This is in j 

t 

addition to the functional difference: between the integrands.) 
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VII. Plots. 


Each plot presents: 

P(«) versus a. 

G(r> versus y. 

H(9) versus 0. 

Plots G~l thru G-6: the range of the independent parameter is 0 to 3. 

Plots H-l thru H-3: the range is o to 10. 

Plots I 1 thru 1 — 3: the range q to 30. 

The distance of the Hydrogen molecule above the paramagnetic 
surface is varied in each series, given as the ratio of z to the 
surface lattice parameter 'a 0 '. 

Each pxot has an 100 subdivisions along the x- or parameter- axis- 
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Abstract . 


The theory for para to ortho hydrogen catalysis on a paramagnetic 
surface, based on the model of a two-dimensional liquid on the surface 
(surface diffusion model) , has been extended to include the effects of 
a surface lifetime *x* . The integral which controls the efficiency 
of para to ortho conversion has been numerically evaluated in terms of 
two dimensionless parameters. The best set of such dimensionless 
parameters was found to be the parameter 'gamma' used in the analysis 
of the case of infinite surface lifetime* and a new parameter: the 
(average) frequency of the para to ortho transition (u ) times the 

Op 

surface lifetime: « . Numerical evaluation of the integral which 
governs the efficiency of the conversion demonstrates some interesting 
features of this model. One of the most significant findings is that 
choosing • 0 p*« ~ an optimal value (about 1.5) leads to a situation 



where the catalysis rate is maximum (for given height above the surface 
of the catalyst) and nearly independent of the parameter (gamma 1 1.5) 
which otherwise controls the efficiency of the rate of catalysis. 
Thus, the residence time 'tau' is of critical importance in achieving 
optimal paramagnetic catalyst performance. 
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I. Method of Analysis. 

The rate of para to ortho catalysis is calculated from the first 
principles result: 


W, 


p — >o 




oO 


u 


d-t 


where 


b 

4. 


M 


H 2 bond length. 

magnetic dipole moment of the proton 
h/2n 


and 


^(t) » a relaxation function 


The general form of ^(t) for any magnetic model for the surface 
catalyst (paramagnetic, ferromagnetic or antiferromagnetic) is as 
follows : 

- a Uj 6D t 

§1-0= 




X x -),» 


U) 


c 0 


1 1 


The quantity (jj)(t) must then be evaluated according to the physical 
model under consideration, and then the indicated integrations and 
summations performed. In practice, the summations are converted to 
integrals, by standard methods used in solid state physics. 
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II. The Two-Dimensional Surface Diffusion Model with Surface 
Lifetime » 'Tau'. 

The surface diffusion model with surface lifetime 'r* is developed 
as indicated in the figure below. 



Surface 



Rate ~ 1 / t 


(Surface Diffusion Layer) 


•z 


(Catalyst Solid Phase) 


The distance 'Z' is assumed to bisect the surface region of average 
distance of the E 2 molecules from the surface plane containing the 
(centers) of regions of space containing the electronic magnetic 
inoniGntib $ in turn locat6d in the surface atoms or molecules . 

According to this model. (j)(t) is determined by evaluating the 
functional forms of the four different terms which control the form of 
the relaxation function. Only one term is different for this model, as 
compared to the model for surface diffusion with infinite lifetime. 
This is the term for the correlation function for movement of the H 2 
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molecule in the surface diffusion layer: 








±/r 


Substituting this relation into the expression for and then 
setting up the expression for the rate of catalysis gives this form of 
the intermediate result for the rate of para to ortho conversion: 


V->o 


f b \ X , Z \ 4 /Ns\ S£S4 i) 

U ) Mf ^ 3 ( w ) 3 


^ 2 . 
lie ll jJrt 



Where : 

D » the diffusion coefficient for a 2 molecular diffusion 
along the surface of the catalyst. 

* ■ the lifetime of the hydrogen molecules j.n the surface 

layer . 

a 0 “ the surface lattice constant 

•‘a " ® a 9 n ©tic moment of the surface atom or molecule 
'A' . 

S - the electronic spin of 'A*. 


i| 

f 


•I 
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III. Analysis in Terms of Dimensionless Parameters. 

The rate of para to ortho catalysis, for the surface diffusional 
model with surface lifetime ’ t ' is given by: 


W. 


p — >o 




where 




and 


4D 

z. 


P 


4-T> T' 

The integrals Gjty.p' and G 2 (r.P> are given by the quantities: 

. <fc* ^ 

f ' t - 


«»<■*> p> = 


asc.s’ 


, 






f 


Gzlv,f0 » 


asc. jP /d 
3 


S - 

a 5 


j 

o 


r ■*. 2, \2. j l 

*p ) * * 


where the coefficient K <~1/Z ) is the same as in the two reports prior 


to this one, and the upper limit of integration is evaluated from: 


q_Z 


«( Z/a^) 




These expressions have the correct limiting form as v — > 00. That 
is. as p — > 0, the result obtained earier for surface diffusion with 
infinite surface lifetime is obtained. 

This is not an optimal set of dimensionless parameters for the 
analysis of the catalytic rate, as numerical calculations of W 

p — > o 

versus * and p quickly show. The problem is that r and p are a bad 
mix of both surface properties and the lifetime on the surface, k much 
better set of dimensionless parameters is this set: 


— * 

CO<»p 

r * * determined by surface properties. 

and 

woptau =* w 0 p*x , determined by the para-- >or tho 

transition frequency and the 
lifetime on the surface. 

This is a more sensible set of parameters, which gives directly a set 
of curves which are readily interpretable. 


IV. Results 


In the figures which follow , we have calculated the para to ortho 
catalysis rate from the expression: 


V->o * * 


^Optau ) ^ 


where the independent parameters are the dimensionless parameters 


_ Z 


and 


42 ) 


woptau » “op** 


The range of woptau is from zero to any large number greater than 1.0, 
for example 0 to 100. We would expect, naively, that the para to ortho 
catalytic rate would peak for a value of <*> ~ unity . 

•he integrals G^(y, woptau) and Gj^y, woptau) are given by: 

V' 2 - , 


G>j [ ^-’op i ^ 




and 


& ^ ^ 
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^ )/ 
* * 




lV> ^o P T\ 






S' - 

$ «■ 




'if 


The range of y is y > 0 . 


IU 


✓ 


^op 1 


*4“ V 



curves . 


The RESULTS are given in the form of two sets of 

S et W; W-l thru W— 9- 

Plots of W p __ >o versus woptau for various values of 'gamma'. in 
general, one can see a maximum in the rate for woptau - 1.5. but the 
results are more informative than this simple rule. 

(1) Starting from the high end of the range of 'gamma', tne catalytic 

rate is found to increase as y decreases, leading towards a maximum 
catalytic rate at r ~ l.j. 

<2) POr r beco ” ln 9 Still smaller, i.e. , < l.s , the entire range of 
catalytic rate versus woptau decreases, i.e. for this condition the 
cataltyic efficiency is gg nexallv reduce . 


Set G: fir-1 thru G-2 . 

In this set, the catalytic rate W 

p — >c 

values of woptau. Conclusions: 


is calculated versus y for cjiven 


(3) For woptau - 1.3 to 1.5 (or some such range), the catalytic rate 
is nearly independent of y for y > about 1.5 


(4) if woptau - 2,4, or 6 
marked decrease versus y. 


(or greater) then the catalytic rate shows a 


V. 


Conclusions on the Optimal Conditions for Para to Ortho 
Hydrogen Catalysis. 

In general this model of surface diffusion which incorporates the 
surface lifetime constant , x* shows: 

1. Distinct optimal properties. 

2 . A distinct l lfl ngg of parameters where the cataltyic rate is 
optimal (i.e for given distance Z above the surface). 

3. A strong dependence of the catalytic rate on the surface life 

time x , expressed through the dimensionless parameter « *x. 

op 

In particular , the optimum lifetime is given by 

x ~ 1 . 5 / <■» 

op 

A typical value of E(ortho) - E(para) is about 50 cm" 1 , from which one 
finds 

“op " 9.4x10 12 radians/sec. 

so that 

■C - 1 .5/(9.4xl0 12 ) - I.6xl0 -13 sec. 

and, generally, : 

t- (optimum) - 10 -13 sec, 
for the optimum surface residence time. 



4. Of course, in all of these observations, the rule that the 
rate - l/Z 6 still holds, so the hydrogen molecule must still come as 
close as possible to the (plane containing the) paramagnetic centers in 
order that the rate coefficient (K) be as large as possible. 
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VII. Figures. 


Set W-l thru W-9. 

Plots versus woptau for given values of the parameter y 
Set G-l thru G-2 . 

Plots versus 'gamma' for various values of 'woptau'. 
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Abstract . 


The theory for para to ortho hydrogen catalysis on a paramagnetic 

surface, based on the model of a two-dimensional gas on the surface has 

been extended to include the effects of a surface lifetime ‘t*. The 

integral which controls the efficiency of para to ortho conversion has 

been numerically evaluated in terms of two dimensionless parameters. 

The best set of such dimensionless parameters was found to be the 

parameter 'alpha' used in the analysis of the case of infinite surface 

lifetime, and a new parameter: the (average) frequency of the para to 

ortho transition ( « ) times the surface lifetime: » . Numerical 

op op 

evaluation of tne integral which governs the efficiency of the 
conversion demonstrates much the same features as for the previous 
model (two dimensional liquid). One of the most significant findings 
is that choosing “ 0 p* T ~ an optimal value (about 1.0) leads to a 


situation where the catalysis rate is maximum (for given height above 
the surface of the catalyst) and independent of the parameter alpha as 
long as alpha is greater than some critical value. We find a(critical) 
on the order of 0 . S . 

Thus, the two fundamental parameters which most control the 

catalytic rate on a paramagnetic surface are: the height. 'Z'. above 

the paramagnetic surface and the surface residence time ’x'. An optimum 

- 13 

catalytic efficiency would exist for x (surface) ~ 10 sec, a 

condition which appears to be difficult to achieve. 

The main conclusions of this study, for the two-dimensional 
gaseous surface layer case, are that the optimal gaseous paramagnetic 
catalyst for the para to ortho hydrogen conversion has: 

1. a minimum distance from the catalytic surface zo the Hydrogen 
molecules in the surface layer, and 

2. the surface lifetime of the H 2 molecules must be as close to 

-13 

10 sec as possible. This is difficult to achieve at a temperature 

-13 

of 120 K, since 10 sec appears to be a distinct lower bound to the 
surface residence time of a hydrogen molecule . 
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I. Method of Analysis. 

The rate of para to ortho catalysis is calculated from the first 
principles result: 


where 


and 




W. 


p — >o 


b 

“l 

K 


£ (i) J $«'«- 


bond length . 

magnetic dipole moment of the proton 
h/2n 




) * a relaxation function 



The general form of ^)(t) for any magnetic model for the surface 
catalyst (paramagnetic. ferromagnetic or antiferromagnetic) is as 
follows : 


5,0 - ® 



Crit)-rw)]\ 

/e 


M- 





U ) (i * 



r v) ' 

/») \ 
\ 4 


The quantity ^(t) must then be evaluated according to the physical 
model under consideration, and then the indicated integrations and 
summations performed. In practice, the summations are converted to 
integrals, by standard methods used in solid 3tate physics. 


II. The Two-Dimensional Gas Model with Surface Lifetime = 'Tau'. 
The two dimensional surface gas model with surface lifetime 'x' is 
developed as indicated in the figure below: 


I 


H 2 Gas 


I 


A 


Rate ~ l/t 


V 



The distance 'Z' is assumed to bisect the surface region of average 
distance of the H 2 molecules from the surface plane containing the 
(centers) of regions of space containing the electronic magnetic 
moments, in turn located in the surface atoms oi molecules. 

According to this model, ^)(t) is determined by evaluating the 
functional forms of the four different terms which control the form of 
the relaxation function. Only one term is different for this model, as 
compared to the previous model for a surface two dimensional gas with 
infinite lifetime. This is the term for the correlation function for 
movement of the H 2 molecule in the two-dimensional surface gas layer: 
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/ C?‘+)-rco^ 

V* 



e 


kT 




-■t/r 


Substituting this relation into the expression for (£(t), and then 
setting up the expression for the rate of catalysis gives this form of 
the intermediate result for the rate of para to ortho conversion: 


V->o 


ibV* / Hi \ 

VK* \~ ) 


S(S-H) 


b_ 


\zz 

) h> Ma ® 




CO 


cos e 


kT 

2m 


\ 


X z 

t 


-*/r 
e d-t 


Where: 

2kT/m 


x 





P 


S 

^c 


mean squared velocity of the H 2 molecule on the surface, 
lifetime of the hydrogen molecules in the surface layer, 
the surface lattice constant 

the magnetic moment of the surface atom or molecule 'A', 
proton magnetic moment, 
the electronic spin of ’A' . 

cut-off value of the lattice wave vector q: n/a . 


J 


III. Analysis in Terms of Dimensionless Parameters. 

This equation for the rate of para to ortho catalysis can only, 
generally, be evaluated by numerical means. We convert this expression 
to a form suitable for numerical evaluation by these transformations: 

1. Replace kT/2m by v„„„ 2 /4. 

rms 


2. Change the integral over 'q' to an integral over x by means of 
the transformation: 


x * 2 q Z 


3. Set the upper limit of the integral over time down from 
infinity to 'T'. where T is chosen to be much larger than the surface 
residence time 'r'. 

4. Change the unit of time in the integral over 't' from seconds 

to the period of the radian frequency: y = w op t ‘ 

5. Introduce the dimensionless constant alpha by the 
substitution: 


* 


w Z / V 

op i,v rms 


Thus we obtain, finally, an expression for the rate of para to ortho 
transition which is: 



and <i> is related to the rotational constant B by 

hw = 2B 

op 

The upper limit to the integral over x is determined as before by the 
relation : 

q„Z = n(Z/a ) 
o 

which requires that we specify the height above the surface in units of 
the surface lattice constant. 

The last step of the transformation of the form of Vf . is to 
c p — >o 

substutute M/2B for 1 /“ 0 p and rearrange the constant terms. The result 



Where the coefficient K has the same meaning as in previous reports: 
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K - 


This expression has the correct limiting form as t --> 00. That is, 
as 1 /t — > 0, the result obtained earlier for the two dimensional gas 
model with infinite surface lifetime is obtained. 

This, now, is an optimal set of dimensionless parameters for the 
analysis of the catalytic rate. This best set of parameters is the 
set : 


and 


U->op 






» determined by properties of the 
2D gas on the surface. 


woptau =* “op* T ' determined by the para — >ortho 

transition frequency and the 
lifetime on the surface. 


This set of dimensionless parameters gives directly a set of curves 
which are readily interpretable. First, however, we must transform the 
double integral above to a form suitable for direct numerical 
integration . 


IV. Numerical Evaluation of the rate of Para — > Ortho Conversion. 



Carrying out a simple rearrangement gives this final form for the 
rate of para to ortho catalysis of a two-dimensional gas on the 
surface, each molecule in the surface having lifetime '*•: 



^ i 

•i 

toopT 

- * -( 

2. 

' * \ 

U* ] 

W K 

s -X 

e 

1 

1 

> 

-P ^ p 

^ i<f ; 


p - >0 48 

f j 
L c 




The numerical integrations are carried out by a subroutine for direct 
numerical double integration. For this purpose, we use the IMSL 
subroutine TWODQ which carries out a two-dimensional iterated integral. 
This subroutine, in turn, approximates the double integral by 

iterated calls to subroutine QDAG which uses a Gauss-Konrad quadrature 

rule for estimation of a single dimensional integral. The over-all 
numerical procedure is slow, but the accuracy can be controlled at the 
expense of slower numerical evaluation of each point on the resulting 
curve of rate versus the parameter value alpha or v/optau. 

Two FORTRAN programs were written for this calculation: 

AM0DEL1 .FOR plots W p __ Q ( in un its of K) versus parameter alpha for 
three chosen values of woptau . 

FM0DEL1 .FOR plots W p — Q ( in units of K) versus the frequency 

parameter woptau for three given values of alpha. 

In each case, the resulting file of relative catalytic rate versus 

alpha (or woptau) was plotted using the TELAGRAF graphics package. 
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V. Results 


In the figures which follow, we have calculated the para to ortho 
catalysis rate from the expression: 


p — > o 


S' -X 

X e. 


u> a p*T 

*■ 

2. 

4 * 

CoS M 

s> r 

We** 


r* 


where the independent parameters are the dimensionless parameters 

tOoD ^ 


>W\S 


woptau = w *r 
op 


The range of woptau is from zero to any large number greater than l.o, 
for example 0 to 100. The upper limit to the 'y '-integral is determined 
by choosing T * 20*x , and the limit to the integral ovar 'x' is chosen 
by setting the height Z equal to one unit of surface spacing, i.e. 


We would expect, naively, that the para to ortho catalytic rate 
would peak for a value of w 0 p*t - unity. This is about what is 
observed . 

The first set of f igures (F-l, F-2, F-3 ) present the results of 

calculating the catalytic rate versus woptau for various values of 
alpha. The results clearly show that the optimum value of alpha 
(for any residence time x) is: 


0.30 


a(optimum) = 

Can this optimum be achieved in practice? To answer this question, we 
calculate the numerical value of o for: 


and 


T = 120°k 


Z * b - 0.74xl0 -8 cm 

which is the covalent bond length of the hydrogen molecule, and 
represents the closest possible distance of the H 2 molecule to a 
catalytic site. Now, at this temperature (in cgs units): 

V rms 2 = 2(8.314X10 7 )(120)/2.016 

so the rms speed of the B 2 molecule on the surface at lOOX would be: 

v rms = "«490 cm/sec 

Nov. using the typical value of %p = 9 . 4x i0 12 radians/sec. we find: 


° - ( 9 . 4x1 0 1 ^ ) ( 0 . 7 4x1 0 8 ) / 9 9 4 90 = 0.70 

Since we expect a larger distance from catalytic site to hydrogen 
molecule than b, this means that an optimum of a (for arbitrary t) 
cannot be achieved, expect for the unlikely condition that the ortho - 
para energy difference were reduced by a factor of 2 - 5 from its 

normal value of about 50 cm '* 1 1 

— s e< »Qnd Sg j; — a£ — f igures — (A- 1.2. 3.4. 5, 6) presents the relative 
catalysis rate versus alpha for various values of wopt.au . Fiqure A-4, 
in particular shows that the maximum catalytic rate is achieved for: 

w op* T “ 1 *° 

which is about what we -naively) expect. However, there is an 
auxiliary condition: That is, if alpha is greater than some critical 
value, then the catalytic rate stays at its maximum efficiency 
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independent of the value of alpha. 


Figure A- 6 pins down this critical value of alpha. This critical 
value of alpha (for optimum tau) is about: 

a(critical) >. 0.5 

fee <*> ** * 1.0 

op 

This condition can be achieved! We note that a typical value of alpha 
(found above) was a = 0.70 (or larger), so this condition could be met 
for a considerable range of distances and ortho-para transition 
frequencies. Therefore, the critical quantity to adjust (besides the 
distance from the catalytic surface... in the l/Z 6 term) is the 
residence time 't*. 
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VI. Conclusions on the Conditions for Optimal Efficiency for Para 
to Ortho Hydrogen Catalysis. 

In general this model of a two dimensional surface gaseous 
layer which incorporates the surface lifetime constant ' t •' shows: 

1. Distinct optimal properties (as for the case of the two- 
dimensional liquid on the surface) . 

2. A distinct range of parameters for which the catalytic rate is 
optimal (for a given distance Z above the catalytic surface). 

3. The optimal conditions can be expressed in terms of the 
lifetime 't* and the dimensionless constant 'a': 


and the condition on alpha is: 


a > 0 , 5 

4. Considerations on the lifetime x . 

1 2 

Now, since we know that a typical value of is 9.4x10 

radians/sec, this tells us that the optimum surface residence time is 
about : 

- 1 3 

x (optimum) - 10 sec. 

While this appears to be a reasonable number, the question is whether 

it is possible to realize so small a value of the surface lifetime. 

To answer this question, we make an estimate of the surface 

lifetime from 3ome elementary mechanics. We assume that the unit of 

time for 'x' is the time for an H 2 molecule to travel one lattice 

constant: a . Wo» take this distance to be two to three Angstroms, or 
o J 
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TV 


a Q - 2.5 A 

Then, at 120 K the rms speed of a hydrogen molecule is 

; v rms ~ 1q5 cm/sec 

and we estimate the unit of time for x by the simple 'time of 

flight' relation: 

a = V T 

o rms 

from which we find: 


t ( unit) 


2 .5x10 8 /io 5 


2.5x10 


-13 


sec 


We expect, moreover, that the value of t on the surface is some 
multiple of this unit of time, for example: 

t = nx(unit) 

where we only know for sure that nil. For example n = 3 to 4 would 
raise the lifetime up to about 10 -12 sec, which would have a 
considerable effect on reducing the catalytic efficiency. Nevertheless, 
we see that r is on the order to the desired magnitude, although 
probably larger . 

5. It is clear that a more thorough study of the surface lifetime 
is in order . since this appears to be one of the most important 
properties of tne catalytic system, the other being the distance Z 
from the surface layer containing the paramagnetic centers. 
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VIII. Figures. 


Figures F-l. Fr2 . F-3 . 

The relative catalytic rate (in units of K) versus «•» *r for a 

op 

set of three values of a. in each figure. 

Figure F-3 establishes an optimum value of a (when considered as a 
function of woptau) , which is : a =» 0.30 

Figures A-i. A-2, a-3. a-4. a-5 . 

The relative catalytic rate ( in units of K) are plotted versus the 
dimensionless parameter a for three values of woptau. in each plot. 

Figure A-4 established that woptau * 1.0 gives the maximum 

relative catalytic rate which is then independent of o for a greater 
than some critical value. 

Figures A-5 and A-6 . 

The critical value of o (when the catalytic rate is considered as 
a function of woptau) is determined by reducing the range of alpha from 
(0,100) to (0,10) and finally (0,2). The peak in figure A-6 (for 
woptau * 1.20) comes at about o * 0.45, This was rounded off in the 
report to o ~ 0.50. 
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THEORY OF PARA TO ORTHO HYDROGEN CATALYSIS. 

March 14, 1990 
R E. Coffman 

Chemistry Department, University of Iowa 


V General Considerations of the Catalytic Rate for Para to Ortho 
Hydrogen Conversion In the Vicinity of the Paramagnetic to 
Ferromagnetic or Anti ferromagnetic Phase Transition. 

Abstract 


Most dense magnetic systems become magnetically ordered at a 
particular critical temperature: The Curie temperature, T c> for 
ferromagnetic systems or the Neel temperature, T^. for antiferro- 
magnetic systems. In the magnetically ordered state, the excitation 
energy for para to ortho hydrogen catalysis may be carried away by the 
magnetic excitations (magnons) which represent long-range oscillations 
of the magnetic spins. This gives rise to two important phenomena: 

1 An increase in the rate of para to ortho hydrogen catalysis in the 
magnetically ordered state, compared to the disordered paramagnetic 
state, and 

2. A critical instability in the rate of catalytic conversion near T c 
or T n with the result that the rate of change of the catalytic rate 
constant (for Increasing Temperature) near the critical temperature may 
become negative. 

Thus a bad choice of catalytic material for the cooling of an exothermic 
process (burning of fuel in a Jet engine) could lead to decreased catalytic 
cooling as an Increase in temperature demanded increased cooling. A 
good understanding of this process and experimental measurements ot 
the catalytic process through an appropriate temperature range would be 
necessary to avoid such a situation 


RATE OF CATALYSIS INCLUDING MANY-SPIN CORRELATION 


The rate of para to ortho hydrogen catalysis was found, in earlier 
reports, from a simple formula which requires calculation of a 
relaxation function and performing at least one Integration: 


where ,s > for an Y magnetic model of the catalytic surface, given by 
the expression: 

h 


X 

For the paramagnetic case, the correlation function of the magnetic 
spins on the catalyst surface which is represented by: 


< S- q (t)Sq(0)> 


has the simple evaluation: 


<S-q(t)Sq(0)> 


I N s S(S*I) 


N N 3 


Where N s is the number of spins, and N is the total number of surface 

atoms of the same type (but not necessarily all paramagnetic). 

Therefore, Ns/N Is the fraction of potential catalytic sites on the 
surface which are paramagnetic. This simple result applies only to the 
paramagnetic case, l.e. the magnetic spins on the catalyst surface have a 
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random character and do not Interact over a long range. 

The correlation function for the magnetic spins on the surface has a 
more complicated form If Interaction with the lattice spin system is to 
Pe included In the rate calculation. Assuming that the spin Hamiltonian 
is pure exchange: 

2< 


(Spins) “^(exchange) 


and then assuming a simplified model (the molecular field orweiss- 
field model) for the exchange Hamiltonian, one finds eventually the 
result (where now N s /N * 1 ): 


<S q (t)S. q (0)> - 

4’VUvt 


4S'WS*> <- 


z 


-..'VUt 


-‘ M ^ 4s v S-> - v cf^ 


where H m is the energy (per spin) of the effective exchange ; eld which 

orders the spins In the ferromagnetic or antiferromagnetic state. 

The result of this changed form of the correlation function is to 
modify the terms In the Integral over time in the basic expression for 
the rate constant giving rise to integral terms of more complicated form: 


oft 

*— 4 ( i. 

< 2 . 


t 



which terms only exist in the ordered magnetic state. (In the disordered 
state, H m - 0). The real part of the expression for this rate will 

therefore have a peak in the vicinity of 

^(o-p) = H m /ft 


which will be modified if other motion or relaxation terms are included 
In the basic expression for the time depe >dence of $(t) 
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Now, the energy of the Interna) field H m is determined in the 
molecular field theory by the Brlllouln function, represented by the 
symbol Bg ; 


3kT c H m S 

H m * BsC ) 

S*1 kT 

which is the formula for the effective field energy, H m ,ln a 

ferromagnetic lattice. Now an antiferromagnetic lattice is made up of 
two sublattices, each of which is ordered by a similar internal magnetic 
field, except that each sublattice Is ordered in the opposite direction to 
Its neighbor. One notes that Is found as a self-consistent quantity, 

and that the value of H m depends on the ratio T /T c . The behavior of H m 
as a function of T/T c is well known (both in theory and experiment!). H m 
Is maximum at T*0 Kelvin, and vanishes at T *T C However, Just below 
the critical temperature <T C or T^) the value of the Internal field, and 

the energy associated with It, Increases dramatically! This therefore 
gives rise to an Increase in the value of the rate constant k(para — > 
ortho) Just below the transition temperature as the catalyst goes over 
from the disordered stata to the ordered state. Physically, this Increase 
corresponds to the additional process of para to ortho excitation energy 
being transferred Into long-range magnetic excitations (magnons) 
which only exist when long-range magnetic order sets in. 


TC 


Temperature 


Figure 1 (Conceptual). Rate constant k( para— > ortho) in the 
vicinity of the magnetic phase transition.. 


One can readily see from the diagram that just below T c or t n , the rate 
of change of the rate constant versus T; 


dk(para — >ortho) 


dT 


a negative quantity 


which can give rise to the difficulties pointed out earlier. 
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NUMERICAL MODELLING OF k(PARA— >QRTHO) IN THE VICINITY OF THE 
MAGNETIC PHASE TRANSITION. 

In the vicinity of the magnetic phase transition, the number and type 
of magnetic excitations, which cooperate with the para— >ortho 
transformation, vary rapidly with temperature. The theory developed to 
this point shows clearly that this process is associated with only one 
quantity (for any model) : the spin correlation function, represented by: 

<S q (t)S-q(0)> 

where the brackets indicate ensemble averages. A means must be 
found to calculate this quantity versus temperature including the effect 
of the phase transition at the critical temperature T c (or T^). If this 
quantity Is then added to the previous expressions for paramagnetic 
catalysis which include the effects of a surface gas or surface two- 
dimensional liquid, then one would irrtve at a theory for the general 
behaviour of k(para— >ortho) In th /iclnlty of the magnetic phase 
transition (as well as away from it!). 

There are two ways to proceed: 

(A) Use the molecular field approximation for the magnetic exchange 
interactloa In this approximation, the exchange Interaction Is replaced 
by an effective (large!) magnetic field interacting with the Zeeman spin 
operator of each magnetic spin. The effective field is represented by the 
Brlllouln function which Incorporates a temperature dependence and by 
solving this function by Implicit means is found to vary with T/T c . This 

is the minimum requirement for a modelling of the average exchange 
interaction over the many spin system, which Incorporates the 
phenomenon of a magnetic phase transition. The concept of magnons, and 
the lattice q' or ’k' vector is missing in this model. This Is the model 
worked out by Ishll. 

(B) Use the magnon concept and work directly with the spin correlation 
function using the Fourier transformed spin operators. This model 
correctly incorporates lattice k-vectors (and magnons!) but does not 
ordinarily Include any average properties of the exchange interaction or 
have dependence on temperature and display a discontinuity in the 
vicinity of the magnetic phase transition. Due to these inconsistencies, 
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this model has not been developed previously since something additional 
Is required In order to make the magnon model work. The extra idea 
arises from the theory of thermal Green's functions which when applied 
to the motion of the magnons and the magnon energy leads to a 
renormalization of the magnon energy In a way which includes a thermal 
average of the z-component of the spin. The renormalized magnons have 
the standard energy at T*0 Kelvin, but have monoton icaliy decreasing 
energy as T-— >T C That Is, they display the correct property of 

vanishing at the critical temperature, and therefore displaying the phase 
transition property. 



Reference: D. Wagner, 'introduction to the Theory of Magnetism", 
Pergamon Press ( 1 972). See Chapter III, Ferromagnetism. 


The idea is this: the conventional magnon energy expression is 
temperature independent, but the renormalized magnons have energy: 

-9(1H - 2s<r<J 0 -J k ) 

where CT * < S z > (which Is proportional to the net magnetization) 
goes too as T — > T c . An equation for <f* is 

4s(f(T) 

cT(T) * i - 

N 

where n^ - the bose-Einstein population of magnons, i.e. 



I 

i 


n k ’ - 

exp^/kT ) - 1 

We can therefore find an Implicit solution for (T(T) at any low 
temperature, defined by T << T c Now, at high temperature, CT” (T) is 
small by definition, so we can find at explicit solution for (T(T) in the 
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vicinity of (Just below) Tc. The solution is easily found to be: 

2kT 

<f (T) - 1 y ^0 ' J k ) 

N * 

where J 0 and J k are exchange coefficients ("exchange integrals"). Since 
the average value of S z , < S z > * (T", * 0 at the critical temperature, 

we can solve for T c : 


KT C N (J 0 - Jj<) 

Thus, the energy of the renormalized magnons > zero as T > T c . 

When incorporated into the theory for the rate constant k.(para > 
ortho) this will give the result that there will be an increase in the rate 
constant just below T c How sharp or diffuse this change may be will 
only be determined by completing the theory for this effect and doing the 
numerical simulations of the rate constant versus temperature. 

The corollary of this Increase In rate constant, will be a negative 
rate of change of k(para — >ortho) versus increasing temperature for T 
Just below T c . However, we cannot yet say Just how far below T c this 

effect will extend. 
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Completion of the Project 


The earliest phase of this project consisted of working through the 
classic calculations of the rate of para to ortho hydrogen catalysis. The 
second stage concerned the working out of a master equation for the rate 
of paramagnetic catalyls, and then working out numerical results for the 
specific models of Hydrogen molecules existing on the surface as a two 
dimensional gas or as a two dimensional liquid. These results, which 
were summarized In four reports (spring, summer and fall of 1989), were 
found to Pe consistent with earlier workers' calculations and gave some 
insight into the role of residence time on the catalytic surface in 
determining the rate of catalysis, as well as the complete set of 
variables which control the net catalytic rate. An important outcome is 
the idea that there is an optimum rate of catalysis for any paramagnetic 
model and a given change in one independent parameter may cause an 
increase in catalytic rate for one catalyst and a decrease for another 
catalyst. (So It is helpful to infer on which side of optimal a given 
catalyst is functioning). 

The unfinished part of this project is to expand the model(s) to 
include the effects of a magnetic phase transition from ferromagnetic to 
paramagnetic, or; anti ferromagnetic to paramagnetic. This is an 
important phenomenon to understand since in the vicinity of this 
transition one may observe negative rates of change of catalyst 
performance with respect to increasing temperature. As with the other 
steps in this study, there is a conceptual problem to be solved first, for 
which I believe I now see the solution. That is, how to Include the 
magnon concept In the modelling In a manner which incorporates the 
phase change behavior in the vicinity of the critical temperature T c 
(ferromagnetic) or T^ (antiferromagnetic). One must then proceed to a 
numerical stage to adequately visualize the significance of these phase 
change phenomena on the para to ortho hydrogen catalyst rate. These 
results should be adequate to explain the general behavior of any 
catalytic system (other than metallic), and give guidance to synthetic 
and other practical studies. 


Theory of Para to Ortho Hydrogen Catalysis 


VI. The Surface Diffusion Model for Paramagnetic Catalysis 
Using Transition State Theory for the Migration of Hydrogen 
Molecules Between Different Catalytic Sites. 

April 12, 1991 

R. E. Coffman 
Chemistry Department 
University of Iowa. 


Abstract 

The theory 'or para to ortho hydrogen catalysis on a paramagnetic 
surface has been developed based on the model of a two-dimensional 
liquid on the surface. This model has now been extended to include the 
effects of a surface lifetime 't' , and a step-size 'd' with an activation 
energy, E d *, using transition state theory for the surface diffusion 

coefficient. In this extended model, the dimensionless parameter 'y 
defined by 


4D 

is interpreted with this formula for the diffusion coefficient: 

D = d^je^’AT 

The important physical parameters which control the catalytic process 
are now: 

z (height above the surface) 
d (diffusion stepsize) and 
a 0 (unit cell length) 
all of which are lengths, and 
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E<j* (activation energy) 

(h/27r)G) 0 p (the first para - ortho splitting) 

both of which are energies, and finally the surface lifetime, Y. 

An analysis of the catalytic rate of conversion of para to ortho 
hydrogen is then developed in terms of four dimensionless parameters 
constructed from the set of dimensioned parameters above. This new set 
of parameters is: 

z/d and z/a 0 

which are length ratios, 

(h/27t)(0 op 

which is a ratio of energies, and finally 

which is a dimensionless product of a frequency times the surface 
lifetime. The integrals for the rate process are then calculated 
numerically and plotted versus the temperature, T ,in degrees Kelvin. 

Careful numerical modelling of this model, for a paramagnetic 
type cataiyst, shows that an optimal choice of parameter values leads to a 
nearly temperature independent catalytic rate over a wide temperature 
range (e.g. 50 to 300 K), but a non-optimal choice of parameters may lead 
to strong temperature dependence over the same range. 

The step size of the diffusion process (through z/d) and the 
relative value of the activation energy of surface diffusion seem to have 
effects primarily related to the shape of the catalytic rate dependence on 
temperature, while the major determinants of the absolute catalytic rate 
are the distance from the surface ( z' through the 'constant' term and the 
z/d dimensionless parameter in the integral) and the surface lifetime. 
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I. Method of Analysis. f 

The rate of para to ortho hydrogen catalysis is calculated from the 


first principles results (derived in earlier reports): 



where 

b = the H2 bond length 

p p » proton magnetic dipole moment. 

h - Planck's constant. 

<t>(t) « a relaxation function 

The general form of <X>(t) for any magnetic model for the surface 
catalysis process (whether the surface :s paramagnetic, ferromagnetic or 
antiferromagnetic) is as follows: 

<D(t) =e-“V £ G 4 (b) X 
I (e-'5{r(0-r(0)]J £ |i<i>(q)p«* (S<j>(t)S«(0)) 

q x,y,z 

The quantity <I>(t) must then be evaluated according to the physical model 
under consideration, and then the indicatea integrations, summations and 
ensemble averages performed. 



II. The Two-dimensional Surface Diffusion Model with 
Surface Lifetime 'Tau'. 

The model for a two dimensional gas on the catalyst surface is 
developed in accordance with the scheme outlined in this figure: 



H2 in Gaseous Phase 



Rate » Js 



Rate - 1 /tau 



Surface Diffusion Layer 


Bulk of Catalyst 


Here, Js represents the rate at wh : ch gaseous molecules col ! :de with the 
surface (calculated from the kinetic theory of gases), and 1/x is the rate 
at which molecules on the surface return to the gaseous phase. The 
distance T represents the height of the surface diffusion layer and is 
defined as the average distance of the hydrogen molecules from the 
surface plane containing the centers of the electronic spin magnetic 
moments. It is these magnetic moments, in turn, which create the 
inhomi aeneous surface magnetic field necessary for flipping the nuclear 
spins of the hydrogen molecules on the surface. 

According to this model for para to ortho hydrogen catalysis, the 
relaxation function <t>(t), is determined by evaluating the functional form 
of the four different terms which de fi r.a the function. The most 
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important term for defining <D(t) in this model is the correlation function 
for the movement of the H2 molecules in the surface diffusion layer: 

( e -iq- [?(t)-r(o)| = e -Dq 2 t e -t/T 

Substituting this relation into the expression for "he relaxation function, 
d>(t), and substitution and evaluation of the other three functions gives, 
finally, the following expression for the rate of para to ortho conversion: 

w - fb'iV 4 _l_ {K\ S(S+1) . x 

^ P ->o-( 2 jm 3 [^2nf ' N ) 3 a 0 2 





q5 e- 2qz dq 


Jo 


COSCOopt 


where: 

D - the diffusion coefficient for the diffusion of hydrogen 
molecules over the catalyst surface. 

x » the lifetime of the hydrogen molecules in the surface layer 
(i.e. the reciprocal of the rate of escape from the surface). 

a 0 - the lattice constant of the crystal face comprising the 
catalyst surface. 

|i A - the magnetic moment of the surface atom or molecule, of 
type ’A’. 

S - the electronic spin of ’A'. 

CDop - the average splitting in energy between the first para and 
ortho pair of states, measured in frequency units : radiaris/sec. 
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III. Para to Ortho Hydrogen Catalysis Rate in Terms of 
Dimensionless Parameters. 

The rate of para to ortho hydrogen catalysis , for this model, may 
be condensed to a simple expression which involves a ’constant' term and 
two integrals which are functions of dimensionless parameters. The 
expression for Wp„ >0 is : 

Wp__>o - K {Gt(y, woptau) + G 2 (y, woptau)} 

where the ’constant 1 K is defined by: 


4tx 2 I 

f N s l 

, S(S+1) | 

■_b_' 

|2 M-p 4 a*" I 

f 1 \ 

Win ) ' 

^ N ) 

1 3 

laoJ 

1 Z 6 1 

I 2 B / 


and the dimensionless parameters which are arguments of the integrals 
G-j and G 2 are defined by: 


y = _2E__ 

r 4D 

and 

woptau = C0 op *T 

The range of the y and woptau parameters is easy to deduce, by studying 
the physical models for the parameters which, in turn, define them. 

The integrals which must be evaluated are defined by: 


r 


oJJl 


Gi(y,© op t) = 


256y 


y7 e - 4 ydy 


y2 + 


J 0 


y \ 2 

j 


+ Y 2 


for the first integral. The other integral is given by: 
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T 


'V 


<*qtZ/2 


Gz(y,co op t) 


25 ef 


y 5 e- 4 Ydy 




These integrals are evaluated numerically on the VAX using the IMSL 
applied mathematics library. 


IV. The Transition State Model for Surface Diffusion. 


In our previous consideration of the integrals which define the 
complex behavior of the para to ortho catalysis process, we considered 
the dimensionless parameter y strictly as a number with a positive range, 
e.g. y > 0. Now, we model y by introducing the transition state theory 
expression for diffusion coefficient *D' expressed in the form: 


D 


-- 1 W 


-Ed*/kT 


where d * the random a alk steplength, kT/h - the frequency factor, and 
E<j* » the activation energy for surface diffusion. When this expression 
for the diffusion coefficient is substituted into the expression for y, and 
some rearrangement is done, the result is: 


K i 

(Z\2 

'(h/27ijoa 0 p 

2 1 

id/ 

L kT J 


e Ed‘/kT 


It is helpfu 1 to define a set of dimensionless ratios, which are 
appropriate to the evaluation of the dimensionless parameters yand 
woptau, and the upper limit of the integral. Clearly, y depends on a ratio 
of lengths (z/d) arid two ratios of energies. The natural unit of the 
activation energy for this problem is clearly the ratio of the activation 
energy to the fundamental para - ortho energy splitting: 


j 



a 


We define this ratio as 'xratio': 




xratio = 


Ed* 


(h/ZjtJoe- 


Then, the expression for the dimensionless parameter y becomes: 



Y 


JX/ZU (h/27tK^ j (h;7l ,v-^ AT 

2 Id/ L kT J 


Now, instead of the two dimensionless parameters that we had before: 


y and woptau 

we now have a new set of dimensionless parameters: 

xratio = (Defined above) 
za ratio - z/a 0 

zdratio - z'd 
woptau - co^/c 

i 

So, the para to ortho catalytic rate will depend on: » 

z and T 

and the four dimensionless parameters above. The sum of the two 
integrals above was therefore calculated and plotted versus the 
Temperature T in degrees Kelvin units. The dependence on 'z' per se would 
require a logarithmic plot, and will be studied separately. It is these 
integrals which define the dependence of the catalytic rate on 
temperature. 
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V. The Results of the Numerical Calculations. 


A. Integral Value versus xratio. 

For this calculation , three of the four dimensionless parameters 
were set to fixed values: 

zaratio - 1 .00 
zdratio - 2.00 
woptau - 1 .4 

(It had been determined earlier that an optimal value of trequency times 
time is wcptau on the order of 1.4). Zaratio = 1.00 means a distance of the 
H2 molecule from the surface equal to the lattice constant, and zdratio = 

2.00 implies a diffusion stepsize of 1/2 of the lattice constant. 

The integrals were then calculated with these fixed parameter 
values, and the ratio of activation energy to para-ortho energy splitting 
was allowed to vary over these ranges: 

xratio - 0.0156, 0.0625 and 0.2500 
xratio - 0.2500, 1.0000 and 4.000 
xratio - 4.000, 8.000 and 16.000 

The quantity plotted is, in each case, the sum of the two integrals 
defining the catalytic rate, versus the temperature T, in degrees K. 

Thus: 


'Integral Value" ■ Gi(y, woptau) + 62(7, woptau) 
The resulting plots are presented on the next three pages. 
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B. Integral Value Versus woptau. 

The previous calculation showed a dependence of catalytic rate on 
the magnitude of the activation energy. In fact, a large activation energy 
(in relation to the para- ortho energy splitting) may be beneficial. A value 
of xratio » 0.25 gave a peak in para to ortho catalytic rate at about 150 K, 
so in the next series xratio was set to this value. 

The three fixed parameters for this series are: 

zaratio « 1.00 
zdratio » 2.00 
xratio =0.25 

and now woptau is chosen as the variable for each plot. This parameter is 
chosen to vary over the ranges: 


woptau « 0.20, 

0.50, 

1.10 

woptau = 1.10, 

1.40, 

1.70 

woptau = 1.70, 

3.00, 

5.00 

woptau =* 5.00, 

10.00, 

20.00 


The resulting plots follow on the next four pages: 
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Integral Value: zaratiosl.OO, zdratio-2 00. xratio=0.25 
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diffmodel2pr.plot 


Integral Value: zaratio=1.00, zdratio=2.00, xratio=0.25 
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C. Integra! Value Versus zaratio. 

The previous calculations and plots show woptau - 1.4 to be an 
optimal value of the transition frequency times residence time, so these 
parameters are assigned fixed values: 

zdratio » 2.00 
woptau -1.4 
xratio - 0.25 

In these calculations the distance from the surface is varied in 
units of the lattice constant, a 0 . (Remember that the 1/z 6 dependence is 
still present in the 'K' term). The results found here measure the effect 
of the hydrogen molecule moving through the field gradient of a para- 
magnetic ion or molecule in the surface layer at various heights from the 
surface. Thus, for z much larger than a 0 , the catalytic rate varies purely 

as 1/z 6 . However, as the molecule get significantly within this distance 
(i.e. z < a 0 , then the rate of catalysis slows down and is slower than 

indicated by the 1/z 6 law). 

The plots shown on the next two pages are for: 

zaratio - 0.500, 1.00, 2.00 

zaratio - 2.000, 3.00, 5.00 




diffmode!3.plot 


Catalytic Integral: zdratio=2.0, woptau=1.4, 
xratio=0.25 
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Diffmodel3P.plot 


CATALYTIC INTEGRAL: idratio=2.0, woptau=1.4, 

xratiosO.25 
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D. Integral Value Versus zdratio. 

The previous set of integral calculations for various values of the 
ratio of height z to the lattice constant a 0 , showed only a small 
dependence on this parameter. (The dependence is accentuated only for 
close distance of approach of the H 2 molecule to the surface). Therefore, 
zaratio is now set to a neutral value: zaratio - 1.00. So, we now fix these 
dimensionless parameters to these values: 

zaratio - 1 .00 
woptau = 1 .4 
xaratio = 0.25 

The integrals were then calculated for the fixed set of parameters 
above and a range of the dimensionless parameter zdratio (the . acio of z/d 
or height of the H 2 above the surface compared to the stepsize of the 


diffusion process). The ranges 

selected 

were: 


zdratio 

= 0.25, 

0.50 , 

1.00 

zdratio 

= 1.00, 

1.50. 

2.00 

zdratio 

* 2.00, 

3.00. 

4.00 


The results are shown on the next three pages: 
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INTEGRAL VALUE: zaratiosl.O, woptau=1.4, xratio=0.25 
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E. integral Value Versus woptau for Optimal Conditions. 

This last set of calculations is for the variation of the integral 
value controlling the shape and temperature dependence of the para to 
ortho catalysis rate assuming that optimal conditions have been chosen 
for: 

zaratio: height z above the surface in units of a 0 . 

zdratio: the distance z compared to the diffusion step size, 
xratio : the relative size of the activa f ion energy for diffusion, 

compared with the para to ortho unit of energy splitting. 

Now, the previous calculation showed that a large stepsize for the 
diffusion mechanism gives a catalytic rate nearly independent of 
temperature. So we optimize the value of zdratio at a value of about 4. We 
also founa in the first part of this integral study (part A), that a ’large’ 
activation energy (compared to the para-ortho energy unit) leads to a 
catalytic rate showing small dependence on the temperature. Thus we fix 
the value of xratio at about 12.0. Therefore, we have selected these fixed 
values of the dimensionless parameters: 

zaratio 1.00 
zdratio * 4.00 
xratio * 12.0 

That is, we have chosen a large value of the activation energy for 
diffusion, since for large values of the activation energy for diffusion 
there is nearly a constant rate of catalysis versus temperature. 

We have then calculated the integral value for these values of 
woptau: 


woptau - 1.10, 1.40, 1.70 

woptau = 1.70, 3.00, 6.00 

The results are shown on the next two pages. 
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Diffmodel2F.plot 


integral Value: zaratio=1.00, zdratio=4.00, xratio=12.0 
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Diffmode!2FF.plot 


Integral Value: zaratio=l.OO, zdratio-4.0, xratio=12.0 
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Diffmodel2G.plot 


Integral Value: zaratio=1.00, zdratios4.00, 
xratso=12.00 
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VI.. Conclusions. 

Magnetic catalysts for the para to ortho conversion may be 
classified in terms of three types of magnetic phases: 

paramagnetic 

ferromagnetic 

antiferrGmagnetic 

For catalysts in the paramagnetic phase, we have derived some 
quantitative results on the behavior of the para to ortho hydrogen 
conversion on the assumption that the hydrogen on the surface behaves as 
a two-dimensional liquid, so that the molecules of H2 migrate from one 
site to another by a Brownian motion mechanism, i.e. by surface diffusion. 

The diffusion coefficient 'D‘ was then modelled using transition 
state theory thus incorporating the stepsize 'd\ and and an activation 
energy for surface diffusion. The numerical expression for the rate of 
catalysis was then analyzed in terms of four dimensionless parameters: 

xratio 

zaratio 

zdratio 

woptau 

and the height from the surface, z, and the absolute temperature, T. The 
first three parameters seem only to affect the shape of the catalytic rate 
versus T, while the most important overall parameters are: 

woptau and z 

Thus, the height above the surface, z, and the residence time, x, seem to be 
the most important determinants of catalytic activity. This is 
particularly evident in the last set of plots given above (E), for which 
the catalytic rate was found to be nearly constant versus the temperature 
for any given value of the surface lifetime. But, increasing the surface 
lifetime was found to measureably degrade the catalyst performance. 


The next step, from here, is to studv a physical model for the 
residence time, x, and the processes that control it. The result will be a 
more or less complete phenomenological theory for this cata.>Jc process. 
One early conclusion: it appears likely that the major differences between 
different paramagnetic catalysts in terms of their catalytic activity are 
the differences in how close the H2 molecules (in the surface layer) 

approach the catalytic site(s), and the length of time that the molecules 
are resident on the surface after once being adsorbed from the gaseous 
phase. Since this (process of adsorption) is a purely physical process, 
with a weak binding energy, such residence times are inherently very 
short. The trick (to obtain optimal catalytic activity) is to shorten the 
residence time to a value such that the parameter designated as 'woptau' 
becomes as close as possible to unity! 
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